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This paper is devoted to study the propagation of Hght beams carrying orbital 
angular momentum in optically anisotropic media. We first review some properties of 
homogeneous anisotropic media, and describe how the paraxial formalism is modified 
in order to proceed with a new approach dealing with a general setting of paraxial 
propagation along uniaxial inhomogeneous media. This approach is suitable for 
describing the space-variant-optical-axis phase plates. 

I. INTRODUCTION 

In the past few years the research field of generation, manipulation and characterization 
of helical beams has attracted considerable interest. Under suitable generation conditions, 
these beams can transport angular momentum {spin and orbital) along their propagation 
direction These beams have triggered a wide range of applications (see the other 

papers in this special issue). Here, we examine the interplay between spin (SAM) and orbital 
angular momentum (0AM) of light upon propagation through anisotropic media, which is 
not only interesting from the fundamental viewpoint, but it also yields novel applications as 
we will expand in the following. 

During paraxial propagation in isotropic homogeneous media, both spin and 0AM are 
conserved quantities |2|, and thereby no exchange of angular momentum between these 
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two is expected. However, in anisotropic media the beam during its propagation can exhibit 
a couphng between SAM and 0AM, even a lossless angular momentum transfer between 
both of them 3[. A proof-of-principle demonstration of spin-controUed-changes in the 0AM 
of circularly-polarized Gaussian beams in the visible domain, using patterned nematic liquid 
crystals, was experimentally achieved in 2006 The strong optical (uniaxial) anisotropic 
properties displayed by nematic liquid crystals were exploited to achieve a controllable space- 
dependent anisotropy of the medium. Moreover, one can reorient the optical axis of these 
media by applying external forces Hence, one can construct optical elements exhibiting 
a space- variant-optical axis, in the sense that the optical axis varies at each point, thus 
allowing the medium to act on polarization differently for each position. A straightforward 
application of this type of media is the creation of non-scalar helical waves, based on spa- 
tially nonuniform polarization transformations, such as the attractive radial polarization or 
azimuthal polarization light beams ds S] and references therein). The other main applica- 
tion is found in the realm of quantum information and communication. Photons can carry 
information encoded in different degrees of freedom: polarization (SAM), spatial modes 
(0AM) and energy. The aim of quantum communication (QC) is to transfer and distribute 
quantum (entangled) states among distant sites (nodes) of networks where they are further 
processed [o], lo|. Devices allowing the unitary (lossless) information transfer from one de- 
gree of freedom to another degree of freedom in the same photon can really improve the 
robustness of the QC network [11]. Anisotropic media are promising tools to transfer the 
entanglement encoded in the polarization of photons to their spatial mode profile and vice 



versa 
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13| . Although previous works have achieved impressive results in this line, it 



is still an open question how to construct truly unitary transformations. Probably more 
complicated space-variant-optical-axis phase plates making use of the anisotropy along the 
radial axis will be necessary to achieve such unitary transformations, in the same line as it 
was anticipated in Ref. 14|. 

So far, due to the interesting applications arisen by the space-variant-optical-axis phase 
plates, it is worthwhile briefly reviewing them. Space-variant-optical-axis phase plates have 
traditionally been fabricated by subwavelength metal stripes jo], and used as polarizers. 
In these, when the period of the grating is smaller than the wavelength of the incident 



light beam, only the zeroth order is a propagating order 



15|, and the gratings behave as 



layers of a uniaxial crystal, whose optical axes are perpendicular and parallel to the grating 
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grooves. This alternative way of creating birefringent components yields polarizing beam 



15|. Exploiting this alternative 



splitters and quarter-wave plates in the visible domain 
approach of fabrication, the Hasman's group has been constructing space-variant-optical- 
axis phase plates to examine polarization manipulations on light. We shall not focus on 
these polarization manipulations in order not to deviate from our main track. For further 
details the reader is referred to Refs. jsl, [l6-18|, where the connection between vectorial 



vortices, scalar vortices, Pancharatnam's phase and adiabatic phase is done. 



for light beams in the midinfrared 
developed space-variant-optical- 



The phase plates of Hasman's group were fabricatec 
spectral domain. However, in 2006 Marrucci et al. 
axis phase plates in the visible domain to prove a unidirectional spin-to-orbital angular 
momentum switching. Spin-to-orbital angular momentum switching has also been observed 



without using space- variant phase plates, for example in the non-paraxial regime |20l . |2l|, 
or a method to control the transfer of spin by an externally applied dc electric field in an 
optically active medium [19[ was also proposed. 

Space-variant-optical-axis phase plates were firstly addressed under the Jones matrix 
formalism jsl, In a previous work, we used a paraxial approach to prove a strong 

modulation in the spin and 0AM changes when varying the traversed optical path [3], but 
the development of such an approach was not presented. Here, we present for the first 
time an alternative approach based on the vectorial paraxial propagation of helical beams 
in space-variant-optical-axis media, which enables us to describe in a remarkable simple way 
the mechanism of spin-to-orbital angular momentum switching. In section [Til we present 



an extension of the theory by Ciattioni et al. [22|, |23| to treat the paraxial propagation 
in homogeneous uniaxial media. In section IIIII we extend the previous formalism to deal 
with the paraxial propagation in inhomogeneous uniaxial media, which constitutes a very 
convenient way to describe the space-variant-optical-axis phase plates. Finally in section 
IIVI we conclude with a particular application of the derived formalism: the spin and 0AM 
variation of Laguerre- Gaussian beams along the so-called g-plate. 



II. PARAXIAL PROPAGATION IN UNIAXIAL MEDIA 



This section is devoted to the propagation of light beams carrying 0AM in optically 
anisotropic media. We will address the interesting case of paraxial propagation along uni- 
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axially homogeneous media, in an analogous way as the previous works of Ciattioni et al. 



22 



24|. 



|23|, but here developing the general case where the optical axis can have any direction 
with respect to the wave propagation. This formalism will be extended to uniaxially inho- 
mogeneous media, presenting simple and compact analytical formulas allowing to deal with 
a wide range of problems. The reader interested in the propagation of cylindrical partial 
waves along bianisotropic media is referred to 

In anisotropic media the Helmholtz wave equation satisfied by the electric field E = 
E(r,t), when considering a monochromatic wave packet E(r, t) = E(r) e"**^*, is: 

^2 

V^E - V(V ■ E) + -^e{uj)E = , (1) 

with oj being the frequency of the monochromatic wave, c the speed of light, and e the relative 
dielectric tensor. Notice that we explicitly include the possibility of a frequency dependence 
for i. In anisotropic media, the third term in Eq. ([1]) couples the field components, whereas 
the second term, which is often neglected in isotropic media, plays a relevant contribution 
here. When e is a symmetric tensor, we can always choose a reference frame where the 
electrical permittivity is expressed in a diagonal form i = diag(e^, ej,, e^,). It is well known 
that two plane-wave families exist satisfying Eq. ([T]) for a given propagation direction, having 
different velocity phase and polarization. These two families are the so-called ordinary and 
extraordinary plane waves. Uniaxial media are materials with the cylindrical symmetry i = 
diag{e±,e±,e\i) around the so-called optical axis (OA), corresponding to e\\, in this case the 
OA is in the z-direction. Defining the wave vector of a plane wave as k = u = ^ n = /cq n 
(u is a unitary vector pointing along the propagation direction, n is the so-called refractive 
index, and n is the vector of refraction), Eq. ([T]) imposes two conditions in the plane wave 
vectors: 

n' = e^, (2) 

!i±^ + ^ = i. (3) 

In uniaxial media is common to define e± = n1 and e\\ = nl, where rio and Ue are the so-called 
ordinary and extraordinary refractive indexes. The dispersion relations ([2]) and ([3]) are the 
main equations to describe uniaxial anisotropic media, which yield two solutions related to 
the refractive index. The first solution, which satisfies Eq. ([2]), does not depend on the 
direction of the wave vector (ordinary waves), at variance with the solution satisfying Eq. 
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([3]), which depends on the wave propagation direction (extraordinary waves). Equation ([3]) 
is vahd in the reference frame where the electrical permittivity is diagonal. In the general 
case, the OA points along an arbitrary direction uqa = (cos sin^, sin ^9 sin 6*, cos^^) = 
cos sin ^ i+sin sin^ j+cos d k, with i, j, k being unitary and orthogonal vectors associated 
with the laboratory frame axes. The dispersion relation ([3]) is modified and can be easily 
obtained considering a rotation of the laboratory frame in order to fix the 2;-axis along the OA 
direction, so that one can resort to Eq. ([3]) with just a simple coordinates transformation. 
Here we have chosen the 9 angle rotation around the y-direction (laboratory frame) and 
then a (p angle rotation around the ^-direction (laboratory frame), obtaining the general 
dispersion relation 

^0 = ^7^ + ir-2 [Kcos^ + kysm^) 

{k^ cos ip + ky simp)'^ (k^ sin Lp — ky cos p>Y , , 

where we have defined 

1 sin^ e cos^ 9 1 _ cos^ 9 sin^ 9 ^ _ ( 1 \ 



nl{9) nl nl ' 'ri'^{9) nl n^ ' An^ \nl nl^ 

When 9 = (the OA is parallel to the 2;-axis) one recovers Eq. ([3]). The two plane wave 
families can be expressed as 

E(±)(r) = Eo(q) exp[t q ■ ± z (k^nl - q^f'^z] , 

EW(r) = Ee(q)exp[2q-rx±2fc,(fc,,yz], (6) 

where q = kx\ + ky] is the transverse wave vector, r_L = xi + ?/j = (r cos 0, r sin 0) is 
the transverse coordinate, the function kzik^^ky) is given by Eq. (|4]), Eo(q) and Ee(q) are 
the ordinary and extraordinary polarization vectors (we will explain below how to calculate 
them). Then, a monochromatic plane wave is completely described by the transverse com- 
ponents kx and ky. Any electric field propagating in an uniaxially anisotropic media, which 
is a solution of the Helmholtz Eq. ([T]), can be written as a superposition of ordinary and 
extraordinary plane waves ([6]), 

E(r) = j d\ {E,(q) exp[? q ■ ± z {klnl - q^'^z] 

+ Ee(q) exp[i q ■ r_L ± i k-^k^, ■ (7) 
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Hence, assuming that a (paraxial) beam is propagating along a well-defined direction, in 
this case ^-direction, this can be perfectly decomposed as in Eq. ([7]), but under the paraxial 
wave approximation {kz ^ kx, ky) we can rewrite Eq. ([71) in a simpler form that describes 
the propagation of any paraxial beam along the uniaxial medium. We start by performing 
the paraxial wave approximation on Eq. (jlj) isolating kz{kx, ky). Therefore we obtain 



kz{kx,ky) — )• —ga; cos 6* sin 6* 



An2 



± koUeiO) 



nm 1 
2nlnlkl 2nlkl 



(8) 



where we have defined 



qx = kx cos ip + ky sin {p , 
Qy = kx sin if — ky cos if . 



(9) 



Now, instead of kx and ky, we use henceforth g^. and qy, which depend of the OA-direction, 
in particular, on the (f variable (which refers to rotations around the z axis). 

Secondly, we need to calculate the polarization vectors Eo(q) and Ee(q) in equation (JT]). 
For the sake of brevity we will only sketch how to derive them. The ordinary polarization is 
perpendicular to the wave vector and the OA, thus Eo(q) oc k x uqa- On the other hand, 
the extraordinary polarization is perpendicular to the Poynting vector and the ordinary 
polarization, Ee(q) oc Eq x S. The direction of the Poynting vector is parallel to the 
gradient of the wave vector surface given by Eq. (j4]). Knowing the direction of the two 
polarizations, we can finally obtain Eo(q) and Ee(q) by imposing the boundary conditions; 
the Fourier transform of the electric field ([7]) at 2; = is 

1 



E(q) 



d^r'^ exp(-z q ■ r'^) E(r'^, 0) 



(10) 



which is equal to E(q) = Eo(q) + Ee(q). Here we provide the results within the paraxial 
approximation, distinguishing two well-defined regimes; when 6 = and 6' 7^ 0. If ^ = 
(the OA is parallel to the 2;-axis), the polarization vectors reduce to 



Eo(q) = — [{klEx - kxkyEy) i + {klEy - kxkyEx) j 
Ee(q) = ^ [{klEx + kxkyEy) i + {klEy + kxkyEx) j 



(11) 

which are in agreement with previous results obtained in Ref. j22|]. On the other hand, 
when 9^0, the polarization vectors become 



Eo(q) = {cos if Ey — sm. if Ex) — sin (y9 i + cos j 
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Ee(q) = (cos V9 i?^ + sin 



^ ; COS sin 

cos 1 + sin 9? J — 



(12) 



An2 

with an explicit dependence on the OA-direction. In this case there is a 2; component in the 
extraordinary part that cannot be neglected even in the paraxial approximation. 

Summarizing, we can expand any electric field into ordinary and extraordinary parts, as 
we did in Eqs. ([7]). Hence, in the paraxial approximation (we still consider the paraxial 
wave propagating along z) the z and xy components of the electric field are 

E^o{r±,z) = e^^°"°^y'Ae*''"""^'i^oE^(q) = e*'''"°'Axo(r^,^), 
^zoir±,z) = , 



X PeE^(q)=e*'"^(^)^A^,(r^,z) 
cos 6 sin ^ 



E.e(rx,^) = e'^°"-(^)-| A(cosv.E.(q)+sinv;4(q)) 



where the projector operators Po and Pe are given below in Eqs. ( !T7j) . Recall that the g^; 
and gy variables, defined in Eq. ([9]), depend on the OA direction. Moreover, we note again 
the presence of a 2; component in the extraordinary field ( |T3i) which cannot be neglected. 
However, in many systems, where the extraordinary refractive index does not differ much 
from the ordinary refractive index, the term is quite small (see Eq. [5]), and we can 

neglect the z component of the electric field (when the OA is collinear or orthogonal to the 
wave propagation, it is a fairly good approximation). Let us examine more features. The 
ordinary part does not depend on the OA direction, as one would expect in uniaxial media. 
However, the extraordinary part has suffered a noticeable change. As the paraxial wave is 
composed by a plane wave times a slowly varying amplitude A_|_e(r_|_, z) (for ordinary waves 
this is A_Lo(r_L, 2;)), we expect that its plane wave-like part will depend on the propagation 
direction, with a refractive index given by ne{6) (for ordinary waves Uo)- The second main 
remark is related to the slowly varying amplitude A_Le(r_L, z), which does not obey the stan- 
dard (isotropic-like) paraxial wave equation, due to the lack of symmetry in the transverse 
plane during the propagation. The slowly varying amplitudes obey 

\ dz ~'~ 2kono^'^) 
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d 

I— + 



dz ' 2kQn1'n? 
where we have defined, 



at. — « sm 6* cos d 



An' 



dr 



A 



(14) 



= X cos (fi + ysiiKp 
Ty = X sin ip — y cos (p 



(15) 



The rotated variables of Eq. f[T5|) simply come from the rotated momenta variables defined 
in Eq. Another interesting point is the occurrence of a new term in the phase of the 
extraordinary wave, proportional to (originated by the walk-off effect). This term is 
responsible for a transversal translation of the extraordinary wave during propagation. But 
most surprising is the fact that projectors Po and Pg are 



h h 



(16) 



when ^ = 0, but in the case when 6* 7^ 0, projectors collapse to 



Po 



sin^ if 



sm ip cos if 



— sm if cos if 

COS^ if 



p. 



COS ip sm if cos if 
sin if cos sin^ ip 



i - iZz{(f)oJz^{-ip) 



(17) 



where TZz{y^) represents a (p angle rotation around the z-axis and o"^ the usual third compo- 
nent Pauli matrix 



cos (p — sm Lf 
sin ip cos ip 

1 



(J. 



-1 



R = 7^^(v?)(3-^7^^(-v?) = 
Notice that projectors ([T7|) are independent of the 9 variable 



cos(2y9) sin(2(y9) 
sm{2(f) — cos(2(/7) 



(18) 
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III. PARAXIAL PROPAGATION IN INHOMOGENEOUS UNIAXIAL MEDIA 

In section [TTl we have examined the propagation of a paraxial wave along anisotropic 
media, considering the OA constant in all the medium, in other words, the medium was 
homogenous and the OA has the same direction at any position of the medium. However, 
what happens if the medium is inhomogenous (i.e., the OA depends on the position)? We 
will solve this problem under certain conditions. 

Before addressing the inhomogenous problem, it is convenient to rewrite Eqs. (fT3|) in 
a suitable form, which will be essential for further calculations. Therefore, recalling the 
Fourier transform (fTOl) . we can substitute it into the transverse electric part of Eqs. (fT3ll . 
obtaining 



where we have used the definitions of the rotated position variables f|T5l) . Notice the simil- 
itude of Eqs. (IT^ with the Fresnel integral in the isotropic case. The Fresnel integral 
describes the evolution of an input paraxial wave, taking into account its initial boundary 
conditions. It is basically composed by the product of the input paraxial wave with the 
Fresnel kernel (propagator). In anisotropic media, integrals ( !T9|) . analogously to the Fresnel 
integral, describe the paraxial wave evolution. Hence, performing the integration with re- 
spect to the variable r', which corresponds to the transverse plane in 2 = 0, we can determine 
the evolution of the wave for any z. First of all, projectors Po and Pe select the ordinary 
and extraordinary part of the input electric field, depending on its polarization. After that, 
each part evolves independently with its corresponding ordinary and extraordinary kernel. 
Notice that the ordinary kernel is invariant under reference system rotations, as the usual 
Fresnel kernel. However, the extraordinary kernel depends on the OA direction, affecting 
thus the wave propagation. We can introduce the following new functions: 




(19) 




(20) 
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which, together with projectors (fT7|) (we are not considering projectors (fT6|) . i.e. the case 
= 0, as the Fourier transform cannot be cast in the same analytical compact expression), 
allow us to simplify expressions ( IT9l) as 



Ex(r^, z) = E^o{r±, z) + E^e(rx, z) = ^f dh'^ 

ZmzJ 



nm 



+ 



7^,(y?)(T,7^,(-^) 



E±(rl). (21) 



Equation ( 1211) has two main terms. The first one is a factor [nl/nl{6)Fe + Fq] times the 
identity matrix, which does not change the input polarization. However, it couples both 
extraordinary and ordinary parts. On the other hand, the second term, [nl/nl{6)Fe — Fo] 
times a rotation matrix, also combines the ordinary and the extraordinary part but, at 
the same time, it varies the input polarization. In the isotropic limit, when Ug tends to 
rio, we find that from [nl/nl{6)Fg + Fo] one recovers the standard Fresnel kernel, whereas 
[nl/nl{6)Fe — Fo] vanishes. 

Equation (ETjl still describes homogenous anisotropic media, with an OA direction given 
by {6, if). However, in an inhomogenous medium, the OA direction can depend on the 
position. Thus, both 9 = 6{r') and if = v'(r') are functions of the position. Inserting 
this dependence into Eq. (12T|) . and remembering that the function Fe also depends on 
the OA direction, we can calculate the evolution of a paraxial wave in uniaxial anisotropic 
inhomogenous media. Of course, we consider that 6{r') and v^(r') are sufficiently smooth 
functions. 

The approach developed here is quite general and allows one to deal with paraxial beams 
propagating along inhomogeneous anisotropic media by using the integral propagator (l2Ti) . 
In particular, the integral propagator fl2Tl) can be used for beams with a certain spatial 
structure carrying 0AM. One can easily apply this integral propagator to describe the 
spin-to-orbital angular momentum switching observed experimentally in Ref. J] by using 
space-variant-optical-axis phase plates. For the sake of completeness, we will address this 
problem with our formalism in the following section. 
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IV. APPLICATION: PROPAGATION OF A LAGUERRE-GAUSSIAN MODE 

ALONG A g-PLATE 

In this section we consider a specific example, the propagation of a Laguerre-Gaussian 
(LG) beam through a space- variant-optical axis medium, whose OA is always perpendicular 
to the propagation direction (this means 0{r') = 7r/2, see the previous section). Here 
we will focus on the, so-called, g-plates, which have the following linear OA dependence 
(p{r') = q (j)' + ao (g and ao are constants). 

Hence, the initial electric field in Eq. ( ET]) is written as E_|_(r', 0', 2; = 0) = (au+ + 6u_) 
LGe^p{r' , (f)'), where Uo- = (i + iai)/\/2 are the circular polarization vectors, and a = ±1 
for right- and left-hand circularly polarization, respectively. The functions LGi^p{r', 0') 
denote the LG modes , which comprisse two terms LG i^p{r' , (j)') = e^^'^' Re^p{r'), where 
Ri^p{r') is a purely radial function depending on the two indices: i and p which account 
for the orbital angular momentum (topological charge) and the radial node number of the 

n fi 

spatial mode, respectively |2[. Using Eq. ( 12T1) we can calculate the propagation of the 
Laguerre-Gaussian beam along the g-plate. The electric field can be separated into two parts 

{r_i,zY'^\ where E^l\r_i,z) preserves the initial polarization 
while that of 'E^l\r_i, z) changes: 

koUo 



E«(r±, = ^Jdh'^{F, + F,){a u+ + 6u_)e^^%p(r') , (22) 
Ef (r^, ^) = ^/ dh'^{F,~F^){b e-''^ ii+ + a e'^^ u.)e'^<^'/2,,,(r') . (23) 



Note that the second term that changes the polarization is also introducing exponential 
terms that include the OA dependence V9(0',ao) = + ao- In this case propagators Fo 
and Fe ( l24l) can be expressed as 



(24) 



where the last approximation in the extraordinary propagator can be done as long as the 
birefringence is not too large, that is — + 11^. The neglected part accounts for a 

small astigmatism in the extraordinary part of the electric field. We will focus on the angular 
part of the integrals ( |22|) and (123!) in order to analyze the spin and angular momentum of 
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the beam along the g-plate. All the integrals in Eqs. ( 122|) and ( 123|) have the form 

'^e^'''>-^i7dr'e^'^(^'+^-''^Re^ (25) 
Aniz Jo ' Jo 



or 



47rzz JO Jo 

where ni and n2 are refractive indexes denoting n^, ne, or {uq + nl)/{2ne)- By resorting to 
the Jacobi- Anger expansion we can easily deal with the angular factor 



^-ihQp.{rr' cos{<f>- (/,')) 



g J, (^rr) e^^(*-<^') . (27) 

k=-oo \ ^ / 



Therefore, by inserting the expansion fl271) in the integrals fl25|) and fl26|) . and integrating 
with respect to the variable 0' (considering the ip{r') dependence in the g-plates), we can 
reduce the integrals to 



or 



^1 k^^^kon,.^^e, |7'rfr'e^^('-^+'-")i?,,,(r')J, (^rA , (28) 



2z 



X 



pdrV^(^^-'^)/?,,(r')J,.,, (^rr') . (29) 



By using integrals fj28|) and fl29l) we can analyze the propagation of the Laguerre-Gaussian 
beam along the g-plate given by Eqs. f l22p and fl23p . The part of the propagation that does 
not modify the polarization of the initial beam is not changing the angular structure of the 
beam neither. Thus, the propagator fl22|) is preserving the polarization and the angular 
structure, in other words, the spin and angular momentum of the paraxial beam remain 



intact. This explains why, in Ref. j4|, there was an experimentally observed part of the 
beam which exhibited unchanged polarization and angular structure. On the other hand, 
the propagator fl25]) does indeed affect the polarization and alters the angular structure too. 
Therefore, a right- or left-handed circularly polarized LG beam changes its polarization to 
left- or right-handed circularly polarized respectively, and its topological charge from i to 
£ + 2q in the case of left-handed circular polarization or to ^ — 2q in the case of right-handed 
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circular polarization. In other words, the variation of spin angular momentum and 0AM 
is —2h and 2qh respectively for an initial right- circularly polarized LG beam, resulting in a 
total variation of angular momentum of 2(g — l)h. For an initial left-circularly polarized LG 
beam the total variation of angular momentum is — 2(g— in complete agreement with the 
experiment performed in Ref. jj]. Moreover we can obtain the radial part of the LG beam 
propagating along the g-plate, which in general will be a superposition of LG modes with 
different radial number p (the radial part of the LG mode is not preserved). Furthermore, 
this formalism is not only able to explain the spin-to-orbital angular momentum switching, 
one can also analytically calculate the angular momentum of light before and after the 
space- variant-optical-axis phase plates, as we first proved in Ref. Sj. 

In this simple example we have showed how resourceful can be the developed formalism 
to address the paraxial propagation in inhomogeneous media. 

V. CONCLUSIONS 

To conclude, we have developed an approach, based on the vectorial paraxial propa- 
gation of helical beams, to deal with homogeneous and inhomogeneous anisotropic media. 
Within this approach, one can calculate the propagation of a paraxial wave through common 

ements such as space- 
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25| . This approach 



anisotropic media, such as polarizers, and more complex optica 
variant phase plates which have shown promising applications [J 
also allows to calculate the angular momentum of light, suitable for describing the 0AM 
changes of the beam during propagation . We have to remark that alternative approaches 



have been developed for some specific cases. Karimi et al. [27j tried to avoid the approxima- 
tion that ^^(r') and ^{v') are very smooth functions in the case of an specific space-variant 
phase plate, the so-called g-plate. However, in order to do so, they neglected the term V ■ E 
in the anisotropic Helmholtz wave equation ([T]). This term is preserved in our approach. In 



a subsequent work, Vaveliuk 28|] showed that the term V ■ E cannot be neglected, providing 
a non-paraxial solution of such complicated anisotropic wave equation. Unfortunately such 
solutions are only valid for g-plates with q = 1. 
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